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The activated surface diffusion of interacting adsorbates is described in terms of the so-called inter- 
acting single adsorbate approximation, which is applied to the diffusion of Na atoms on Cu(OOl) for 
coverages up to 20% in quasielastic He atom scattering experiments. This approximation essentially 
consists of solving the standard Langevin equation with two noise sources and frictions: a Gaussian 
white noise accounting for the friction with the substrate, and a white shot noise characterized by 
a coUisional friction simulating the adsorbate-adsorbate collisions. The broadenings undergone by 
the quasielastic peak are found to be in very good agreement with the experimental data reported 
at two surface temperatures 200 and 300 K. 



PACS numbers: 68.43.Jk, 05.10.Gg, 68.35.Fx 



I. INTRODUCTION 

The diffusion of atoms, molecules, or small clusters on 
metal surfaces plays an important role in many techno- 
logical and industrial applications. Molecular beam epi- 
taxy, heterogeneous catalysis, or the fabrication of nanos- 
tructures are examples strongly affected by the kinetics of 
diffusion. Among the different experimental techniques 
utilized to analyze the diffusion and adsorbate dynamics, 
quasielastic helium-atom scattering (QHAS) is a gentle 
and inert technique commonly applied to fast diffusioufi 
QHAS has been applied to different metal/metal and ad- 
sorbate/metal systems, the diffusion of Na adatoms (at 
different coverages) on Cu(OOl) being one of the most ex- 
tensively systems studied so far ,^i^iii^i^iL^i^ i i o 1 1 1 1 1 ^ At low 
coverages, this system has been theoretically described 
within the so-called single adsorbate approximation and 
analyzed in terms of the motional narrowing efject,'^''^ 
and also within the Kramers' turnover theory^ and the 
dephasing theory.^ In particular, the motional narrowing 
effect governs the broadening of the lineshapes as a func- 
tion of the friction, the parallel momentum transfer, the 
lattice structure, and the jump dynamics. Generalized 
expressions for the lineshape have been recently given 
under these different conditions 

When dealing with higher coverages, absorbate- 
adsorbate interactions can no longer be neglected. 
Then, pairwise interaction potentials are usually intro- 
duced into two-dimensional Langevin molecular dynam- 
ics (LMD) simulations^ where the number of coupled 
equations to be solved is usually very high (it increases 
as 2N, where TV is the number of adparticles considered, 
typically of the order of 400-500). On the one hand, first, 
this translates into a high computational cost. Second, 
the numerical results issued from LMD simulations are 
difficult to interpret since analytical treatments are not 
easily implementable. On the other hand, from a physi- 



cal viewpoint, the broadening (measured as the full width 
at half maximum) of both the quasielastic ( Q) peak rul- 
ing the diffusion process and the T-mode peaks related 
to the low frequency motions of the adsorbate {frustrated 
translational modes), is not well reproduced when com- 
pared with the experimental datai^ Recently, it has been 
shown§. that a good agreement with the experiment can 
be achieved using simple models where the adsorbate is 
allowed to also move perpendicular to the surface. 

Aimed to provide a theoretical and numerical alterna- 
tive to the standard procedure at intermediate coverages, 
we have recently proposed the so-called interacting single 
adsorbate (ISA) approximation. ^'^'^^'^^ Within this ap- 
proach, diffusion is described by using only one standard 
Langevin equation, which is characterized by the follow- 
ing three contributions: (1) the deterministic, adiabatic 
potential V, which models the adsorbate-substrate inter- 
action at T = 0; (2) a Gaussian white noise, Rcit), ac- 
counting for the lattice vibrational effects that the surface 
temperature induces on the adsorbate; and (3) a white 
shot noise, which stands for the adsorbate-adsorbate col- 
lisions and replaces the pairwise (dipole-dipole) interac- 
tion potential generally considered. In this way, a typical 
LMD simulation involving N adsorbates is substituted 
by the dynamics of a single adsorbate, with the action 
of the remaining 7V-1 adsorbates being described by the 
random force associated with the white shot noise. This 
noise describes the adsorbate-adsorbate interaction by a 
series of random pulses within a Markovian regime (i.e., 
pulses of relatively short duration in comparison with 
the system relaxation). This interaction is, therefore, de- 
scribed in terms of a coUisional friction which depends on 
the surface coverage. With this simple stochastic model, 
a better agreement with the experimental data for cov- 
erages up to around 20%, approximately, is obtained. 
Although further investigation at microscopic level and 
calculations from first principles are needed, this simple 
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stochastic model at moderate coverages is also able to 
provide a complementary view of diffusion and low fre- 
quency vibrational motions, described by the peaks at 
or around zero energy transfers (very long time dynam- 
ical processes), respectively. This could be understood 
because any trace of the true interaction potential seems 
to be wiped out due to the relatively large number of 
collisions taking place at very long times. 

The organization of this paper is as follows. In Sec.HIl 
we introduce the stochastic model used in the ISA ap- 
proximation as well as its connection with the observable 
magnitudes in QHAS experiments. In Sec. IIIH the work- 
ing model considered for the diffusion of Na on Cu(OOl) 
is described, and numerical results are presented and dis- 
cussed. Final remarks and conclusions are summarized 
in Sec. lIVl 



II. INTERACTING SINGLE ADSORBATE 
APPROXIMATION 

In analogy to scattering of slow neutrons by crystals 
and liquidsfi^ii^ the observable magnitude in QHAS ex- 
periments is the differential reflection coefficient, 



d^7^(AK,^) 
dQduj 



(1) 



This coefficient gives the probability that the He atoms 
scattered from the adsorbates on the surface reach a cer- 
tain solid angle Q with an energy exchange hiu — Ef — Ei 
and wave vector transfer parallel to the surface AK = 
— Ki. In Eq. ([T]), rid is the concentration of adpar- 
ticles; J- is the atomic form factor, which depends on 
the interaction potential between the probe atoms in the 
beam and the adparticles on the surface; and S{A'K,uj) 
is the dynamic structure factor or scattering law which 
gives the Q and T peaks and it provides a complete in- 
formation about the dynamics and structure of the ad- 
sorbates through particle distribution functions. Exper- 
imental information about long distance correlations is 
obtained from the scattering law when considering small 
values of AK, while information on long time correlations 
is provided at small energy transfers, huj. 

The dynamic structure factor in Eq. ^ can also be 
expressed as 



S{AK,uj) = ^ I e-"^* I{AK,t)dt, 
2tt J 



(2) 



where 



/(AK t) EE (e"*'^^-[^(*)"^(°)l) = (^(,-i^i<fo-"^Kit')dt''^ 

(3) 

is the so-called intermediate scattering function. In this 
function, the brackets denote the ensemble average over 
trajectories R(i) running on the surface, and uak is the 
velocity of the adparticle projected onto the direction 
of the parallel momentum transfer, AK {AK = ||AK||). 



Following the standard formulation used in neutron scat- 
tering theory, this function can be split up into two 
contributions: Is {self function), which corresponds to 
the average of trajectories of the same adparticle at two 
different times, and Id {distinct function), which cor- 
responds to the average of trajectories of two different 
adparticles at different times. The Fourier transforms 
corresponding to / and Ig give what are called the co- 
herent and incoherent scattering laws, S'(AK,ti;) and 
^^(AK,^), respectively. 

In QHAS experiments, and with interacting adsor- 
bates, coherent scattering is always obtained. The corre- 
sponding theoretical interpretation of this type of scat- 
tering is usually carried out in terms of Vineyard's con- 
volution approximation,-!^ where the distinct pair corre- 
lation function is expressed as a convolution of the self 
pair correlation function. This approximation is known 
to fail at small distances, where the surface lattice be- 
comes important. Whereas in neutron scattering many 
attempts to improve the convolution approximation have 
been developed, within the QHAS context very little ef- 
fort has been devoted to achieve this goal. At finite cov- 
erages, one usually distinguishes between two diffusion 
coefficients)^ the tracer diffusion constant {Dt) and the 
collective diffusion constant {Dc). Dt refers to the self- 
diffusion process and focuses on the motion of a single 
adsorbate. On the contrary, Dc is related to the collec- 
tive motion of all adsorbates which is governed by Fick's 
law. In either case, a Kubo-Green formula relates Dt or 
Dc with the velocity autocorrelation function of a single 
adsorbate or with the corresponding for the velocity of 
the center of mass, respectively. 

In the ISA approximation, the distinction between self 
and distinct time-dependent functions does not apply, 
and Eqs. ([!]) and ([2]) still hold. If the so-called Gaussian 
approximatioiilS' is invoked, some analytical treatment is 
possible. This theoretical approach leads to a very simple 
and easy manner to interpret the underlying dynamics 
issued from the numerical Langevin simulations. Under 
this approximation, the intermediate scattering function 
can be expressed as a second-order cumulant expansion 
in AK as 



I{AK,t) w e- 



AK^ J*(t-t')Cit'}dt' 



(4) 



In Eq. ([4]), C{t) is the velocity autocorrelation function, 
defined as 



C(r) ^ {v{0)v{t)) 



(5) 



along AK. The Gaussian approximation is exact when 
the velocity correlations at more than two different times 
are negligible, thus allowing to replace the average acting 
over the exponential function by an average acting over 
its argument. 

The motion of an adsorbate under the action of a bath 
consisting of other adsorbates on a static two-dimensional 
general surface potential can be well described, in the 
Markovian approximatioufi by the standard Langevin 
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equation 

R(i) = -vMt) + FMt)] + S^s{t) + SKait), (6) 

where R = {x,y) represents the surface lattice points and 
F — — Vy is the deterministic force per mass unit, with 
V{x, y) being the surface interaction potential with pe- 
riods a and h along the x and y directions, respectively. 
The substrate excitations lead to the random force per 
mass unit 'R.Q{t) on the adatom which has features of 
a Gaussian white noise. The shot noise Rs(i) is given 
by a series of impacts, and we have assumed sudden 
adsorbate-adsorbate colhsions (i.e., strong but elastic) 
and after-collision effects relax exponentially. Moreover, 
the probability for collisions follows a Poisson distribu- 
tion. No correlation has been assumed between both 
types of noise as well as between the noise sources as- 
sociated with each degree of freedom. 

In Eq. rj — J + X, where 7 is the frictional damp- 
ing coefficient resulting from the nonadiabatic coupling 
to the electronic and vibrational excitations of the sub- 
strate, and A is the number of collisions per time unit 
or collisional friction. By using the elementary kinetic 
theory of transport in gases^^ and the Chapman-Enskog 
theory for hard spheres, a simple relation can be foundii 
between the collisional friction coefficient A and the cov- 
erage at a temperature T, 



A = 



6pe ksT 



(7) 



where fcs is the Botzmann constant, a is the unit cell 
length of an assumed square surface lattice cell and p is 
the effective radius of an adparticle of mass m. 

Finally, it is worth stressing the fact that, at long 
times, there will not be free diffusion because the effects 
of the two stochastic forces are dominant. This is the 
diffusive regime, where mean square displacements are 
linear with time. 



2kBT 
mq 



t = 2Dt, 



(8) 



D being the diffusion constant. In general, we can express 
Eq. (0) as 



i? = ^lhn l(|R(t)-R(0)p}, 



(9) 



which is known as Einstein's law, the value D — ksT /mrj 
for the diffusion coefficient being the Einstein relation. 
Note from Eq. ([8]) that D increases by lowering the total 
friction rj and by increasing the surface temperature. 



III. RESULTS AND DISCUSSION 

The diffusion of Na atoms on the Cu(OOl) surface can 
be considered as a prototype in QHAS since a lot of ex- 
perimental and theoretical work can be found in the lit- 
erature. The coverage = 1 corresponds to one Na 



atom per Cu(OOl) surface atom- or, equivalently, a = 
1.53 X 10^^ atom/cm^; a — 2.557 A is the unit cell length 
and p — 2 K has been used for the atomic radius of Na. 
The surface friction we have considered in our simulations 
is taken from Ref. [3 to be 7 = 0.1 ujq = 2.2049 x 10"^ 
where is the harmonic frequency associated with the 
periodic surface potential. Frequencies and times are 
here given in atomic units. We are going to analyze four 
coverages for two surface temperatures: 200 and 300 K. 
On the one hand, ^Na = 0.028 and 0.064, where the single 
adsorbate regime holds; on the other hand, ^Na ~ 0.106 
and 0.18, where the interaction between adsorbates can- 
not be ignored. These coverage values correspond to 
A = 3.34 X 10-^ (?7 = 2.53 x lO^^), A = 7.64 x 10^^ 
(r? = 2.98 X 10^5), A = 1.26 x 10"^ (r? = 3.46 x 10"^) and 
A = 2.15 X 10^^ (77 = 4.68 x 10"^), respectively. 

For the adsorbate-substrate interaction potential V , in 
this work we have considered the non-separable model 
proposed by Ellis et al.,^ whose form is 

V{x, y) = Vo{x, y) + [x, y) -f V2{x, y). (10) 

In this expression, the first term is a simple separable 
cosine potential, 

Vq{x, y) = [2 - cos(27ra;/a) - cos(27ry/a)] , (11) 

with a the lattice constant of the Cu(OOl) surface and 
Vb =41.4 meV; the second term, which reads as 

Vi{x,y) = -Ay^e"''^t^/''"(™+^/^)l'+['^/""("+^/^^l'^ 



m.n 



(12) 

with A — 2Vo and b = 11.8, is added to produce a low- 
ering of the potential barrier at on-top sites according to 
the observations; finally, the third term is also a nonsep- 
arable part which serves to alter the curvature near the 
minima and vary the difference between the potential at 
the minima and the bridge positions, 

V2{x,y) = 7r2CVb^[(a;/a-m)2 + (y/a-n)2] 

m,n 

X exp[-{x/a - m)"^ ~ {y/a ~ n)"^], (13) 

with C = —0.2. Two scattering azimuths, [100] and [110], 
diagonal and parallel to the reciprocal space lattice in the 
first Brillouin zone, respectively, are considered here. 

In the ISA approximation, Einstein's law, given by 
Eq. ([9]), is satisfied and, therefore, the diffusion constant 
D decreases with the coverage. As said above, the origin 
of the friction arises from the collisions among adsor- 
bates and the influence of the substrate phonons. This 
behavior is clearly observed in Fig. [TJ where two types 
of calculations of the diffusion constant are plotted as a 
function of the total friction (and, therefore, of the cover- 
age, since the surface friction coefficient is kept fixed in all 
our calculations) and for a surface temperature of 200 K 
along the azimuth [100]. The full circles correspond to 
the values obtained from Eq. ^ after performing the cor- 
responding Langevin numerical simulation. The dashed 
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FIG. 1: Numerical (circles) and theoretical (dashed line) val- 
ues for the diffusion coefficient _D as a function of the total 
friction at T = 200 K along the azimuth [100]. 



line represents Einstein's relation with a fitting constant 
instead of the actual value fcsT/m; the value obtained for 
such a fitting constant would correspond to an effective 
mass 1.007 times the Na mass. As seen, the agreement 
is fairly good. 

The broadening observed in the Q peak with increasing 
coverage can be attributedi^iiiii^ to the respective fast 
decay undergone by the intermediate scattering function. 
The maximum value of the experimental peak widths is 
found to be increased by a factor of 3 as the coverage is 
increased from 0.028 to 0.18. In Figs. [Ua)-(d), the re- 
sults of our model for the Q-peak width have been plotted 
and compared with the experimental ones for = 0.028, 
0.064, 0.106 and 0.18 at T = 200 K [black squares (sohd 
lines)] and T — 300 K [red squares (dotted hues)] to- 
gether with the experimental ones (black and red trian- 
gles, respectively) along the azimuth [100]. To compare 
with, in each case a few LMD results are also shown (blue 
stars). As can be seen, LMD simulations reproduce the 
general trend, but predict a smaller increase in the broad- 
ening of the Q peak.-^ However, the agreement of the ISA 
results with the experiment is fairly good up to = 0.106. 
This broadening has to be understood in terms of what 
we have termed as hound trajectories^ associated with 
the motion of the adsorbates inside the potential wells: 
when adsorbates display this type of motion, the com- 
plex exponentials in Eq. ([3]) keep the coherence necessary 
for the intermediate scattering function to decay much 
slower than when the adsorbates arc moving on the sur- 
face {running trajectories) . This effect is directly related 
to the increment of r], which is equivalent to consider an 
increment of the coUisional friction A (which increases 
with 0), since we have chosen 7 as a constant substrate 
friction in all our numerical simulations. The same trend 
could be observed in the lineshapes associated with the 
T-mode peaksJ^^iiii^ 

As mentioned above, the agreement of our model with 
the experiment is at about 9 = 0.106. It is clear then that 
this value should be regarded as a limit for our model. 
Note that with high coverages, the adsorbate-adsorbate 
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FIG. 2: (Color online.) Numerical (squares) and experimental 
(triangles) dependence of T on l^K at T = 200 K [black (solid) 
line] and T — 300 K [red (dotted) line] along the azimuth 
[100]. Different values of the coverage are considered: (a) 
e = 0.028, (b) e = 0.064, (c) 9 = 0.106, and (d) 9 = 0.18. 
To compare with, in each case a few LMD results (blue stars) 
are also shown. 
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FIG. 3: Numerical (squares) and experimental (triangles) de- 
pendence of r on AK at T = 200 K along the azimuth [lIO]. 
Two values of the coverage are considered: (a) 6 = 0.028 and 
(b) e = 0.18. 



interaction starts to play a more relevant role in the dif- 
fusion dynamics. The average motion of the adsorbate 
becomes slower and it feels the force exerted by its neigh- 
bors for longer times, which could end up in a sort of 
collective motion involving a number of adsorbates. This 
would give rise to a breakdown of the Markovian ap- 
proximation and, therefore, taking into account memory 
effects. 

Broadening with the coverage can also be observed 
along the [110] azimuth. As far as we know, no LMD 
simulations have been performed in this direction. Our 
numerical results (black squares) at two coverages {9 — 
0.028 and 0.18) and the experimental data (triangles) at 
9 = 0.028 are displayed in Figs.[3la) andlS^b). In partic- 



ular, the results for the highest value of the coverage are 
considered to be predictive. 



IV. CONCLUSIONS 

A full but simple stochastic description of the line- 
shapes for activated surface diffusion and low-frequency 
vibrational motion for adsorbate/substrate systems with 
increasing coverages has been presented and discussed. 
A simple Langevin equation within the ISA approxima- 
tion is used to calculate the adparticle dynamics, where 
a Gaussian white noise accounting for the effect of the 
adsorbate-substrate friction and a white shot noise sim- 
ulating the friction due to the collisions among adsor- 
bates, with no correlation between them, are considered 
to describe microscopic interactions. In our opinion, 
the present numerical results for the broadening of the 
Q peak, which are obtained with no fitting parameter, 
are in better agreement than previous LMD simulations, 
mainly due to the role of the coUisional friction intro- 
duced in our stochastic model. We are well aware that the 
barrier diffusion should be modified with the coveragci^ 
However, our goal here was only to study the collisional 
friction effects. Obviously, this type of studies cannot 
replace further investigations at microscopic level, and 
calculations from first principles of this dynamics need 
to be done. The complementary view we propose in this 
work should remain valid because the lineshapes of the 
diffusion and low-vibrational frequency motion analyzed 
at very long times, where the number of collisions is rel- 
atively large (or where the statistical limit is applicable), 
seem to wipe out any trace of microscopic interactions. 
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